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Preface 


This document is an update to JPL Document 80-93 dated February 15, 1981. Its two purposes 
are to present an overview of the accomplishments of the Deep Space Network in the field of 
radio science research, and to provide to the space science community information on existing as 
well as currently planned funded Network capabilities. As part of the baseline capability, we have 
included the implementation being carried out for the Cassini mission to Saturn and Titan. It is 
planned that future radio science experiments and the anticipated requirements associated with 
them will be addressed in a separate document by the same authors. 

This document has had the benefit of review by members of the radio science community. 
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Abstract 


Radio science experiments use radio links between spacecraft and sensor instrumentation that is 
implemented in the Deep Space Network. The deep space communication complexes along with 
the telecommunications subsystem on board the spacecraft constitute the major elements of the 
radio science instrumentation. Investigators examine small changes in the phase and/or amplitude 
of the radio signal propagating from a spacecraft to study the atmospheric and ionospheric 
structure of planets and satellites, planetary gravitational fields, shapes, masses, planetary rings, 
ephemerides of planets, solar corona, magnetic fields, cometary comae, and such aspects of the 
theory of general relativity as gravitational waves and gravitational redshift. 
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I. Introduction 


The Deep Space Network is a telecommunications facility managed by the Jet Propulsion 
Laboratory of the California Institute of Technology for NASA. Its primary function is to provide 
two-way communications between the Earth and the spacecraft exploring our solar system. To 
carry out this function, it is instrumented with high-power transmitters, low-noise amplifiers and 
receivers, and appropriate monitoring and control capabilities. It consists of three complexes 
situated at approximately equally spaced longitudinal intervals around the globe at Goldstone, 
California; Robledo near Madrid, Spain; and Tidbinbilla near Canberra, Australia. Two of the 
complexes are situated in the northern hemisphere while the third is situated in the southern 
hemisphere. Waff (1993) provides a historical retrospective on the Network's development into 
the present configuration. 

The Network comprises four subnets, each of which contains one antenna at each complex. The 
four subnets, illustrated in figure 1, are the 70-m diameter subnet (where the stations are 
designated DSS 14, DSS 43, and DSS 63), the 34-m diameter standard subnet (DSS 12, DSS 42, 
and DSS 61), the 34-m diameter high-efficiency subnet (DSS 15, DSS 45, and DSS 65), and the 
26-m diameter subnet (DSS 16, DSS 46, and DSS 66). Another subnet is currently in the 
construction phase. It will consist of three 34-m beamwave guide antennas which, in addition to 
X-band (8.4 GHz) electronics, will also provide efficient communication at Ka-band (32 GHz). It 
may also transmit at that band to support the Cassini mission radio science experiments. 

These deep space communication complexes along with the telecommunications subsystem on 
board the spacecraft constitute the major elements of the radio science instrumentation. Radio 
science investigators examine the small changes in the phase and/or amplitude of the radio signal 
propagating from a spacecraft to an Earth receiving station in order to study the atmospheric and 
ionospheric structure of planets and satellites, planetary gravitational fields, shapes, masses, 
planetary rings, ephemerides of planets, solar plasma, magnetic fields, cometary comae, and 
aspects of the theory of general relativity like gravitational waves and gravitational redshift. 
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Figure 1. The Deep Space Network 1992 Configuration 




II. The Network as a Radio Science Instrument 


In addition to communicating with the deep space missions, the Network generates accurate radio 
science data observables. The following is an overview description of the various subsystems at 
the deep space stations that affect the radio science data. Figures 2 and 3 are block diagrams of 
these subsystems and illustrate how they interact with each other as well as with a spacecraft radio 
instrument. Figure 2 shows the example of the Mars Observer mission and figure 3 shows the 
example of the Galileo mission, without the spacecraft but with additional details of the ground 
data system. Table 1 lists selected characteristics of the Network's stations applicable to radio 
science. 

The large parabolic surface of the antenna focuses the incoming RF energy onto a subreflector. 
The subreflector is adjusted in position to optimize the transfer of energy to the rest of the system 
housed inside the antenna and elsewhere in the complex. In transmission, RF energy is focused on 
the subreflector which is adjusted to optimize the transfer of energy to the primary reflector and 
from there to the target spacecraft. 

There are several possible ways to keep the massive antennas properly pointed at the spacecraft. 
One method is "CONSCAN" where the acquired signal is conically scanned by the antenna. The 
feedback from the closed-loop receiver is provided to an algorithm which compares the scan 
pattern to the received signal level and compensates to point the scan center at the apparent 
direction of the spacecraft signal. 

In cases of excessive signal dynamics or low received signal levels, as in many occultation 
experiments, CONSCAN cannot be used and instead the antenna is "blind pointed." Pointing 
angle predictions are computed for each antenna using an ephemeris provided by the spacecraft 
navigators. Before the tracking period begins, the predicts are transformed into coordinates for 
the antenna. An antenna-mounted computer adjusts the coordinates for subreflector position and 
refraction, and adds a systematic error correction. 

With the 70-m and 34-m standard subnets, either of two receiver channels can receive S-band (2.3 
GHz) or X-band (8.4 GHz) frequencies. Currently these stations can transmit only at S-band. One 
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of two circular polarizations is permitted for each frequency band. With the 34-m high-efficiency 
subnet, only one receiver channel is available and is optimized for X-band frequencies on 
downlink and uplink (although S-band reception is also permitted). The 34-m antennas have a 
low-power transmitter (20 kW) available while the 70-m stations also have the high-power 
transmitter (400 kW). 

Once the signal has been focused by the front-end of the antenna, it is channeled to a low-noise 
amplifier (typically a maser) and then to a receiver. There are two general forms of reception in 
the Network. One is closed-loop reception which provides, via feedback loop in the circuitry, the 
capability for rapid acquisition of a spacecraft signal and telemetry lockup. In order to achieve 
acquisition within the shortest time, the closed-loop receivers may also be predict driven. The 
carrier signal and ranging spectra are passed to equipment which counts the Doppler cycles and 
computes the spacecraft range and then transmits the data to JPL in real time. The Doppler and 
range data types, as well as signal strength information, are processed by the tracking system 
centered around a metric data assembly and a sequential ranging assembly. 

The metric data assembly controls all the tracking system functions at the station. It counts the 
Doppler cycles and then formats and transmits them to JPL. The ranging assembly measures the 
round-trip light time of a radio signal traveling from the ground station to the spacecraft and back. 
As shown in figure 4, a coded signal is modulated onto the carrier and transmitted to the 
spacecraft where it is detected and transponded back to the station. The received code is then 
correlated with the transmitted code and the round-trip light-time delay of the signal is measured 
(see Kinman, 1992, for a review of Doppler tracking of planetary spacecraft). A range observable, 
r k , is a measure of round-trip light-time (RTLT) as recorded by the station clock (see Asmar and 
Kursinski, 1991): 


r k = ST(t k ) - ST(t k - RTLT) 

where ST(t k ) is the time as kept by the station clock at true time t k . A Doppler observable, d k , is 
defined in terms of differenced range and the Doppler count time: 


dk “ ( r k-l ~ r k) I (Jk+1 ” l k) 
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The other type of signal reception is called open-loop. The open-loop receivers do not have a 
feedback loop with which to track the received signal; instead they rely on frequency predicts to 
remain tuned to the incoming signal. The signal is downconverted from RF to IF in the antenna 
and then transmitted to the complex's signal processing center where it is further downconverted 
to VF. The RF-IF downconverters in the 70-m stations have four IF channels that may be 
recorded simultaneously. They are the right and left hand circular polarizations of the S- and X- 
band signals (i.e., XRCP, SRCP, XLCP, and SLCP). The RF-IF downconverters in the 34-m 
stations have two IF channels: S-band and X-band of either polarization. 

A radio science processor digitizes and records the output from the IF-VF downconverters. The 
data as well as the supporting header information are transmitted to JPL in real time or recorded 
on magnetic tape at the station. The processor tunes the open-loop receiver based on predicts and 
allows for frequency and time offsets. It can also snap tune between the coherent and non- 
coherent frequency channels transmitted by the spacecraft. The radio science instrumentation at 
the Network also transmits a spectral signal indicator display to the Network Operations Control 
Center and the Mission Support Area. The displays include the received signal spectra at the 
station, which are used in real time to validate the radio science system data and monitor the 
progress of the experiment. 

The predicts used for open-loop receiver tuning, closed-loop signal acquisition, as well as antenna 
pointing are generated by the Network based on the latest trajectory information available for the 
spacecraft. Models of the planetary atmospheres are also included in the predict-generation 
software for most of the planets in the solar system. For occultation experiments, it is critical to 
the process of tuning the open-loop receivers to use predicts with the most accurate atmospheric 
models and the latest trajectory information available from the mission navigators. 

All of the functions necessary to operate and control a complex are centrally located in the 
monitor and control subsystem. Through this subsystem, the complex receives, stores, and 
distributes the information needed to operate the stations (e.g., schedules, pointing and tuning 
predicts). Status information from each subsystem is routed through the monitor and control 
subsystem to the staff on site and at JPL. Event alarm messages as well as operator logs are 
processed and displayed through the monitor and control subsystem, where an operator can 
monitor and control the entire group of equipment (called a link) required to support a given 
tracking pass. 
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Several subsystems at the stations require frequency and timing reference. The typical reference is 
a hydrogen maser standard; cesium standards act as backup. These highly stable clocks are critical 
to the precision of the radio science data. Current standards provide stability near 10"^ at 1000 
seconds integration time. 

In addition to the radio science data, the Network also makes available ancillary data on the media 
calibration in the Earth atmosphere. Included are ionospheric calibration data obtained using 
Earth-orbiting satellites and tropospheric calibration data obtained using surface measurements of 
pressure, temperature, and relative humidity. Water vapor radiometers may be used in the future. 

The Ground Communications Facility supports all the communication requirements of the Deep 
Space Network including the radio science system. This facility is distributed among different 
physical locations at JPL and at the complexes. Data lines are used to transmit schedules, 
sequences of events, and predicts to the complexes from the Network Operations Control Center 
before a pass, and system status and configuration information from the complex to the operations 
center during a pass. If necessary, data can also be played back over these data lines at the end of 
a pass. In real time, the computers of either the radio science support team or the Advanced 
Multimission Operations System at JPL capture and display all radio science data relevant to the 
experiments in progress. Voice nets keep the entire operation team in communication by linking 
the operators at the complexes with the controllers at JPL as well as the project and radio science 
personnel. 

The Radio Science Mission Support Area is the control center for radio science experiments. It is 
staffed by the flight projects' multi-mission radio science support team for the purposes of 
conducting operations and monitoring the progress of the experiments. Requirements for 
configuration for a specific radio science experiment are provided by the radio science support 
team, representing the flight project and the principal investigators. The configuration typically 
includes selecting the appropriate stations, open-loop equipment filter bandwidth and sampling 
rates, antenna pointing method, Doppler sampling rate, and ranging parameters, based on a 
detailed analysis unique to that experiment (see, for example, Kursinski and Asmar, 1991). 

Detailed description of the radio science instrumentation is available in several documents. 
Interested investigators and engineers are advised, as a starting point, to refer to the latest 
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revisions of the following JPL documents: Deep Space Network/Flight Project Interface Design 
Handbook (1988), Deep Space Network/Systems Requirements Detailed Interface Design (1991), 
and Asmar and Herrera (1993). 
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Figure 2. Example of end-to-end radio science instrumentation block diagram 
for the Mars Observer spacecraft and Deep Space Station 














10 


Figure 3. Example of Galileo radio science ground instrumentation block diagram 
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Figure 4a. Two-way coherent Doppler measurement 
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Figure 4b. Phase-locked loop at Deep Space Station 
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Figure 4c. Range measurement 

Figure 4. Doppler and Range Measurements (from Kinman, 1992) 

PLL: Phase-locked loop TWTA: Travelling wave tube amplifier 
LNA: Low noise amplifier VCO: Voltage-controlled oscillator 
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III. Radio Propagation and Occultation Investigations 


Investigations of planetary atmospheres, ionospheres, rings, and magnetic fields using radio 
propagation techniques have been conducted by almost every planetary mission and are planned 
for many future ones. In these techniques, radio frequency transmission from a spacecraft 
occulted by a planet (or one of its satellites) and received on Earth probe the extended 
atmosphere of a planet. As a result, the radio link is perturbed in phase and amplitude and this 
perturbation is measured on the ground and converted into an appropriate refractivity profile 
depending on the object of the investigation. From the refractivity profile, information is derived 
about the electron distribution in the ionosphere, temperature-pressure profile in the neutral 
atmosphere, or particle size distribution of the ring material surrounding a planet, in the case of a 
ring occultation. The refractivity profile is obtained only after the Doppler effect due to the 
geometrical flight path has been computed and the observed data have been corrected for it. The 
application of this technique to solar occultations is described in a different section. 

The history of radio propagation studies goes back to the early 1960s. Von Eshleman of Stanford 
University first proposed the method in 1962. Independently, a JPL team led by Dan Cain, 
conducting a study of the effect of refraction in the Earth's atmosphere on the accuracy of the 
counted Doppler, realized the possibility of applying the sensitivity of Doppler phase 
measurements to the study of the atmospheres and ionospheres of other planets (Kliore et al., 
1964). The theoretical technique was developed largely by Fjeldbo (1964). The occultation 
techniques were demonstrated for the first time in 1965 when Mariner IV flew by Mars; they were 
used to determine the salient features of the Martian atmosphere, showing that the atmosphere is 
predominantly CO2 and that the surface pressure is less than one percent of that of Earth, or an 
order of magnitude less than had been previously believed (Kliore et al., 1965; Fjeldbo and 
Eshleman, 1968). 

A descriptive overview of radio propagation experiments in the outer solar system with the 
Voyager spacecraft is provided by Eshleman et al. (1977) and Tyler (1987). An overview of Mars 
Observer radio science experiments is provided by Tyler et al. (1992). The Galileo experiments 
are described by Howard et al. (1992), and the Ulysses Io plasma torus occultation experiment by 
Bird et al. (1992c). Mariner, Viking, Pioneer, and Magellan experiments are cited below. 
Yakovlev (1985) describes radio propagation experiments on Soviet missions. 
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After the first experiment was performed by the Mariner IV mission to Mars, which provided the 
first measurement of the Martian surface pressure and temperature as well as a measurement of 
electron density in the Martian ionosphere, Venus was the next planet probed via radio 
occultation during the 1967 Mariner V mission (Kliore et al., 1967; Fjeldbo et al., 1971). This 
experiment also revealed interesting new information about the temperature and pressure profiles 
in the neutral atmosphere and electron density distribution in the ionosphere. The results of this 
experiment led to the realization that the Soviet Venera 4 spacecraft did not reach the surface of 
Venus (Kliore and Cain, 1968; Eshleman et al., 1968). 

In order to eliminate the loss of data due to closed-loop receiver lock-up time at the occultation 
egress and to provide higher stability instrumentation, the Deep Space Network augmented the 
receiving equipment at the deep space stations with open-loop receivers (see Section II). The 
new equipment was called the radio science system and has been used by all occultation 
experiments since the system's implementation. 

During the late 1960s and early 1970s, the inner planets continued to be studied by radio 
scientists. Mariners 6 and 7 flew by Mars (Kliore et al., 1969), then the Mariner 9 orbiter made 
repeated measurements of Mars (Kliore et al., 1972, 1973); Mariner 10 visited Venus and 
Mercury (Howard et al., 1974a,b; Fjeldbo et al., 1976). The two Viking Orbiters and Landers 
provided large amounts of important occultation data during a mission to Mars which began in 
1976 and ended in 1982 (Fjeldbo et al., 1977; Lindal et al., 1979). 

From December 1979 to October 1992, the Pioneer Venus spacecraft orbited Venus and provided 
important data on the structure of the ionosphere (Kliore et al., 1979b,c) and atmosphere (Kliore 
and Patel, 1980; Cimino et al., 1980; Woo et al., 1980), and their changes in response to the solar 
cycle (Kliore et al., 1991; Kliore and Luhmann, 1991). 

In 1973, the Pioneer 10 spacecraft passed by Jupiter, allowing the first occultation measurement 
of an outer planet. Analysis of these data disclosed the existence of an atmosphere of Io (Kliore 
and Woiceshyn, 1976; Kliore et al., 1975). In 1974, Jupiter was probed by the Pioneer 11 
spacecraft. An occultation experiment was performed during the Pioneer 1 1 encounter of Saturn, 
and this experiment yielded the first measurements of the atmosphere and ionosphere of Saturn 
(Kliore et al., 1980a,b). 
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The Voyager 1 and 2 spacecraft flew by Jupiter in 1979 and Saturn in 1980 and 1981. Voyager 2 
also flew by Uranus in 1986 and Neptune in 1989. Some of the findings of the Voyager radio 
science experiments are summarized below (Tyler, 1987; Tyler et al., 1989). During the Voyager 
occultations of the atmospheres of the four giant planets, pressure level information was obtained 
at about 1.0, 1.3, 2.3, and 6.5 bars at Jupiter, Saturn, Uranus, and Neptune, respectively (Lindal, 
1992; Tyler et al., 1989). Information can be obtained from these temperature-pressure profiles 
about the near-adiabatic lapse rates below the tropopause, at about 100 mbar pressure, as well as 
the considerable structure in the stratosphere. These profiles also provide information on the 
atmospheric constituents, specifically the ratio of hydrogen to helium. Figure 5 shows the 
occultation geometries and figure 6 shows the temperature-pressure profiles for the four outer 
planets encountered by Voyagers 1 and 2. On Neptune, Saturn, and Jupiter, the depth of the 
measurements was limited by microwave absorption. In the case of Uranus, the depth was limited 
by the spacecraft trajectory which did not allow link bending angles greater then 3° to be 
observed. The measurements at Neptune were conducted near 60° north latitude, and the Uranian 
atmosphere was probed near the equator. The data from Saturn and Jupiter were obtained at 
latitudes of 36° north and 12° south respectively. The Jovian temperature distribution may have 
been affected by vertically propagating atmospheric waves (Lindal, 1992). 

Voyager also studied the ionospheres of Jupiter, Saturn, and Triton. The ionosphere of Saturn is 
complex and highly variable with location and time of day. Similar variations with height and 
location were also reported for Jupiter and appear, from preliminary analysis, to be present on 
Uranus as well (Tyler, 1987). 

Eshleman et al. (1979b) describe the Io torus experiment at Jupiter. Levy et al. (1981) studied the 
dispersive Doppler measurement of the electron content of the Io torus from the Voyager 1 radio 
occultation data. Using the occultation data, studies of planetary shapes as well as the magnetic 
fields of Jupiter and Saturn (Hinson, 1984) were also possible. 

The radius of Titan was determined to be 2575 ± 0.5 km during Voyager l's near-diametric 
occultation. It was also shown that nitrogen is the major constituent of Titan's atmosphere. The 
temperature-pressure profile is shown in Lindal et al., 1983; it terminates at a surface temperature 
of 94 ± 0.7 K. Scintillations were observed at all altitudes below about 100 km, reflecting 
surprising dynamics in the atmosphere. 
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The atmospheric radio occultation techniques were extended to planetary ring occultation 
experiments. During the Voyager 1 encounter of Saturn, the ring occultation experiment led to 
observations spanning the entire width of the classical ring system. The results are discussed in 
Tyler (1987); also see publications by Marouf, Rosen, and Gresh in the references. Marouf et al. 
(1982) describe the theory of ring occultation experiments. 

During the Ulysses flyby of Jupiter in February 1992, a path which provided a gravity assist to a 
polar orbit of the Sun, an occultation experiment was conducted on the Io plasma torus. The 
spacecraft transmitted signals at 2.3 and 8.4 GHz. The signals were received, recorded, 
differenced, and integrated in order to derive the columnar electron density of the Io plasma torus 
(Bird et al., 1992c). The measurements agreed with contemporary models of the torus based on 
Voyager data, implying that the amount of gas injected from Io is similar to that observed during 
the Voyager era. On the other hand, the torus seemed to be less extended out in the centrifugal 
equator, implying a smaller plasma temperature than predicted. 

In 1991 and 1992, the Magellan spacecraft also conducted occultation experiments to study the 
atmosphere of Venus, obtaining highly accurate profiles of atmospheric refractivity and 
absorptivity (StefFes et al., 1993; Jenkins et al., 1993). These profiles, in turn, provided profiles of 
ionospheric density, atmospheric temperature, pressure, and sulfuric-acid-vapor abundance. A 
spacecraft limb-tracking maneuver, plus the high effective isotropic radiated power, made it 
possible to probe the atmosphere to depths of 35 km at the 3.6 cm wavelength and to 33 km at 
the 13 cm wavelength. 
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Figure 5a. Voyager occultation geometries at Jupiter, 
Saturn and Uranus. Figure gives plane-of-sky view of 
Voyager occultation experiments. Earth is shown to 
scale (from Tyler, 1987) 



Figure 5b. Mew from Earth of Voyager occultation of 
Neptune (from Lindal, 1992) 


Figure 5. Voyager occultation geometries 
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Figure 6. Temperature profiles for the giant planets derived from radio 
occultation data acquired with the Voyager spacecraft (from Linda), 1992) 






IV. Solar Corona and Solar Wind Research 


An overview of radio sounding techniques for studies of the solar corona has been recently 
presented by Bird and Edenhofer (1990). The first use of the Deep Space Network in solar 
coronal research consisted of a spectral broadening experiment performed by Goldstein (1967) as 
the Mariner 4 signal passed within 0.6 deg of the solar disk. Goldstein and colleagues used the not 
yet completed Goldstone 64-m station to receive the very weak Mariner 4 signals (weak because 
the spacecraft antenna was not pointed at Earth) via an open-loop receiver driven by an oscillator 
programmed with the predicted spacecraft received frequency profile. The 100 kW uplink was 
provided by the Goldstone 26-m antenna. In 1968, Goldstein (1969) duplicated the experiment 
with the superior conjunction of Pioneer 6 and obtained results similar to those from the Mariner 
4 experiment. The radially dependent spectral broadening obtained by Goldstein and colleagues 
reflected the now better understood relationship between spectral broadening and integrated 
coronal electron density. 

More recently, Woo (1978) obtained extensive near-Sun spectral broadening measurements of the 
Helios 1 and 2 spacecraft, and used these data in an examination of spectral broadening. A 
substantial set of Helios and Pioneer spectral broadening data has been empirically compared to 
an electron density (N e ) model of the following form, where r is the radial distance and A and B 
are fit constants: 


N e = At- 6 + Br~ 2 3 


For spacecraft near solar conjunction, range measurements as a function of solar offset distance 
provide an excellent method of mapping electron density in the solar corona, provided all other 
competing effects, such as trajectory errors, unmodeled spacecraft forces, and relativistic effects 
can be properly accounted for. Muhleman et al . (1971) first used this method during the solar 
conjunctions of Mariners 6 and 7. Using ranging data acquired at the Goldstone 64-m station, 
Muhleman and colleagues (Muhleman et al. , 1977) were able to model the coronal electron 
density as 


N e = At- 6 + Br~ 8 
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where e = 2.05 for Mariner 6, and 2.08 for Mariner 7, a functional form which is in good 
agreement with earlier eclipse measurements. A similar experiment was performed some years 
later using Helios, which produced results consistent with those of Mariners 6 and 7. Bird and 
Edenhofer (1990) summarized in figure 7 the work of various investigations on the solar wind 
velocity as a function of solar distance. 

The effect of columnar charged particles on the group velocity of electromagnetic waves is given 
by Koehler (1968) as 

r = Kf^ 2 ! N e dR 

where 

r = observed range delay, m 
K = 40.3 

f = carrier frequency, Hz 
N e = electron density, m -3 
R = signal path, m 

Another radio sounding experiment was performed in November 1968, when the Pioneer 6 
spacecraft was occulted by the Sun. Pioneer 6 transmitted a linearly polarized carrier, which 
allowed the first experiment to be performed using a human-made signal source to determine the 
Faraday rotation of the solar corona (Stelzried, 1970). The quasi-longitudinal approximation for 
Faraday rotation is given by (Stelzried et al., 1970): 

O = Qf- 2 J NeB L dR 

where 

Q = Faraday rotation, deg 
f = signal carrier frequency, Hz 
Q = 1.3548 x 10 6 
R = signal path, m 
N e = electron density, m~ 3 

Bl = longitudinal component of solar magnetic field, T 
Measurement of coronal Faraday rotation, therefore, provides information about both electron 
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density and the solar magnetic field (Rusch and Stelzried, 1972). To support this experiment, the 
Goldstone 64-m station was equipped with a remotely rotatable, microwave linear feed (Reid et 
al., 1973). A closed-loop polarimeter (Ohlson et al., 1974) was implemented and used to 
automatically track the orientation of the received signal polarization. This system maximized the 
received signal strength and yielded precise signal polarization data. Observations during 
November-December 1968 yielded information on both solar coronal transient events (Levy et al., 
1969) and the background solar corona. Research using the polarization tracking capability at the 
Network's three 64-m stations (now 70-m stations) was continued with the Pioneer and, 
especially, the Helios spacecraft (Volland et al., 1977; Bird et al., 1977). The Helios Faraday 
rotation data were also used in studies of the magnetic field in coronal mass ejections (Bird et al., 
1985) and the mean background coronal field during solar minimum conditions in 1975-76 
(Patzold et al., 1987). Dennison et al. (1978) made further use of these results with a 
measurement of the gravitational deflection of a polarized signal propagating through the solar 
corona. Another Faraday rotation experiment was conducted with the Magellan spacecraft in 
1992. Figure 8 from Howard et al. (1992) shows coronal Faraday rotation data recorded 
simultaneously at widely separated Deep Space Stations. 

Using single frequency (S-band) Doppler data acquired during the 1975 solar conjunctions of 
Pioneer 10, Pioneer 11, and Helios 1, Berman and Wackley (1976) demonstrated that the 
observed Doppler fluctuations provide an excellent measure of solar wind turbulence. This was 
significant in that Doppler noise was already being automatically recorded for purposes of ground 
tracking system performance monitoring, and hence represented a free source of scientific data for 
coronal studies. Using this data type, it was demonstrated (Berman et al., 1981) that the mean 
level of Doppler fluctuations is proportional to integrated electron density. 

A far more powerful tool for measuring electron density is available when dual-frequency carriers 
are transmitted from a spacecraft. In this case, differenced range measurements made at the two 
frequencies provide a direct measurement of columnar electron density that is not contaminated 
with the multitude of other factors previously mentioned. The Viking spacecraft were equipped 
with this dual-frequency capability, and the 1976-77 and 1978-79 Viking solar conjunctions 
provided excellent electron density measurements (Tyler et al., 1977; Muhleman and Anderson, 
1981). Dual-frequency ranging observations acquired during solar occultations of the Voyager 
spacecraft in 1985 and 1988 were reported by Anderson et al. (1987) and Krisher et al. (1991) 
respectively. Kinman and Asmar (1988) used the 1985 Voyager 2 solar conjunction to study the 
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two-way coherent Doppler error due to the solar corona inducing phase scintillations on the 
uplink. 

Utilizing dual-frequency (S-band minus X-band) Doppler data from the Mariner 10 spacecraft. 
Woo et al. (1977) demonstrated that the columnar phase fluctuation spectral index was 
approximately -2.6: 

W 0 (v) <x v -2 - 6 

where 

W 0 = columnar phase spectral density, rad 2 /Hz 
v = fluctuation frequency, Hz 

This work made a significant contribution to the gradual resolution of the controversy about the 
value of the spectral index, which originated in the early days of (natural source) interplanetary 
scintillation observations. More recently, Woo and Armstrong (1979) have studied the variation 
of the spectral index with radial distance using Viking solar conjunction data and have concluded 
that a significant change occurs in the spectral index at a solar offset distance of approximately 20 
solar radii. 

In addition to the radial dependence of Doppler noise, the highly dynamic solar wind causes 
temporal variations in Doppler noise. Statistical studies show that the occurrence of Doppler 
noise transients depends strongly on the phase of the solar cycle. While the rate of the transients 
is about one every four days during solar maximum conditions, it drops by a factor of 3 to 4 
during solar minimum (Woo, 1988; Armstrong et al., 1992). Comparison of Doppler noise with in 
situ plasma measurements indicates that most of the transients are interplanetary shocks during 
solar maximum (Woo and Schwenn, 1991), while during solar minimum, some of the transients 
are corotating high speed streams (Woo and Armstrong, 1992). Examples of this comparison are 
reproduced in figure 9, showing the time histories of Helios 1 in situ measurements of proton 
temperature T p , proton density n p , proton mass flux density n p v p , and Pioneer Venus rms Doppler 
scintillation a D during April 23-28, 1981. The three Doppler transients T2, T3, and T4 detected 
by Pioneer Venus at 15 to 18 solar radii were observed later at 170 to 174 solar radii by Helios 1 
as shocks S2, S3, and S4, respectively. 

The intensity of the Doppler noise transients, as characterized by the enhancement above a 
background Doppler noise level, is proportional to, and hence indicative of, the speed of the 
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interplanetary shock. The largest transients seem to represent the fastest shocks, which decelerate 
strongly while propagating away from the Sun, and are associated with flares (Woo and 
Armstrong, 1981; Woo et al., 1985; Woo, 1988), On the other hand, the smaller transients seem 
to represent shocks with relatively low speeds. These move uniformly through the interplanetary 
medium and are more often associated with the eruption and ejection of cool filamentary material 
from coronal quiescent regions (Woo et al., 1990). 
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Figure 7. Summary of various investigations of the solar wind velocity as a function 
of solar distances (from Bird and Edenhofer, 1990) 
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Figure 8. Faraday rotation data recorded simultaneously at widely separated Deep Space Stations 
(from Howard et al., 1992) 
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APRIL 1981 


Figure 9. Time histories of Helios 1 in situ data (top 4 panels) and Pioneer Venus 
Doppler scintillation data (bottom panel) showing same features; see text 
(from Woo and Schwenn, 1991) 


V. Celestial Mechanics Investigations 


Perhaps the earliest radio science application of Deep Space Network capabilities occurred in the 
area of celestial mechanics. Previous to the era of deep space exploration, calculation of the 
physical properties of the planets (mass, size, shape, and orbit) had to be conducted via remote 
optical astronomical techniques. The ability to make accurate distance and velocity measurements 
of the deep space probes that were perturbed by a planet's gravitational field allowed in situ 
measurements of planetary constants for the first time. Gravitational properties of the planet's 
satellites can also be determined. In the large planet/satellite systems like Jupiter and Uranus, for 
example, satellites can create motions of the planet which then produce a detectable effect on the 
spacecraft. 

A fundamental property of a planet is its gravitational field, which encompasses many aspects of a 
planet's composition and motion. This field is represented by an nth degree and order (D&O) 
two-dimensional spherical harmonic series. The planetary mass scales the harmonics of the field. 
The ratio of the Sun's mass to the mass of a planetary system is an important celestial mechanics 
parameter, used in calculating the ephemerides of planetary bodies. Additional lower degree terms 
in a field describe the dynamic oblateness and the rotational axis, as well as place constants on the 
deep interior structure of a body. Higher degree terms describe in increasing detail shallow and/or 
localized gravity features of the planet's surface, including mountains and craters. 

Knowledge of a gravitational field, along with a determination of the surface topography and the 
planet's volume (achieved through optical imaging, radio science occultations, or onboard 
spacecraft altimeters), produces the mean density of a planet. The information about the internal 
structure of the planet, when combined with an optical or radar map of surface features, can place 
constraints on models describing the differentiation of interior layers and the chemical 
composition and physical state of the interior. In turn, these models can be used in arguments for 
different theories of solar system evolution. For example, different ideas on satellite origins or 
planetary condensation can be strengthened or ruled out based on the abundance of certain 
materials in moons or asteroids. 

The use of radio science in celestial mechanics experiments began with the 1962 Mariner 2 
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mission to Venus and the 1962 Ranger mission to the Moon. Using the Network-generated 
Doppler navigation data from the Mariner mission, Anderson and Warner (1966) were able to 
make order-of-magnitude improvements over the previous, earthbound, determinations of the 
masses of the Moon and Venus. With the Ranger mission, estimates were made of the masses of 
the Moon and Earth and an estimate was made of the offset between the lunar center of gravity 
and center of figure, Subsequently, Mariner 4 (Mars), Mariner 5 (Venus), and Mariner 9 (Mars) 
provided the mass and first estimates of low order gravitational harmonics for Venus and Mars, 
respectively (see Anderson, 1974, for a review), while analysis of Doppler measurements from the 
Lunar Orbiter and Apollo missions led to the discovery of large positive gravity anomalies 
("mascons" for mass concentration) on the Moon (Muller and Sjogren, 1968). 

Mercury was also observed during this period, by two Mariner 10 flybys (1974 and 1975). 
Analysis of the data yielded the mass and dynamical oblateness of Mercury (Esposito et al., 1978; 
Howard et al., 1974a). A recent re-analysis of the same data has produced a mass accurate to 
0.4%, second order coefficients of the gravity field, and identification of a gravity anomaly in the 
vicinity of the spacecraft's closest approach point (Anderson et al., 1986a). 

Venus has received intense scrutiny in the time since the first Mariners. The data volume from the 
long-term orbiting spacecraft was equivalent to thousands of single flybys. Most notably, the 
Pioneer Venus Orbiter, in operation since 1979, and the Magellan radar mapper, in orbit since 
1990, have produced a wealth of information, making the Venus gravity field the best known 
planetary field after the Earth's, and many studies have been made using the data from this period. 
The Pioneer Venus Orbiter and Magellan data have been combined in at least one study, 
producing a 60th degree and order spherical harmonic gravity model and 120th degree and order 
spherical harmonic topography model, which are the most accurate to date (Konopliv et al., 
1993). Previous studies have produced estimates of fields to 50th degree and order, using only 
Pioneer Venus Orbiter data (Nerem et al., 1993). 

Various gravity anomalies have also been identified on Venus. These features of the gravity field 
are especially interesting because the topography of Venus has been mapped with such high 
accuracy. Correlating the gravity and topography data sets provides invaluable information on the 
dynamics of the Venusian crust. Areas having received particular attention include the Ishtar 
Terra plateau (Sjogren et al., 1984), the Beta and Atla Regions (Sjogren et al., 1983), and 
Aphrodite Terra (Black et al., 1988), with many other regions also explored (for examples see 
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Smreker and Phillips, 1991). 


Changing orbital characteristics of these spacecraft, including periapsis, latitude, and eccentricity, 
allow new features of the planet and its gravity field to be identified and allow different 
perspectives on old information. Magellan has suspended radar mapping and, in September 1992, 
lowered its periapsis altitude to 182 km for a Venus gravity campaign, which is expected to last 
until May 1993. At that time, aerobraking will result in circularization of the orbit for additional 
gravity coverage. 

Mars has also had the benefit of orbiting spacecraft, with the start of the Mariner 9 orbit in 1971 
and the two Viking orbiters arriving in 1976. The best gravity model is a 50th degree and order 
field (Nerem et al., 1993); its best resolution is near the equator due to the large eccentricities of 
the orbiters and periapsis latitude being nearly equatorial. So the poles, for example, are hardly 
characterized at all. As with Venus, Mars data have exhibited many interesting gravity anomalies, 
whose corresponding surface features have been recorded by the Viking cameras. These features 
include Olympus Mons, which is the largest known gravity anomaly in the solar system (see 
figures 10 and 11), the mascons in the Isidis basin (Sjogren, 1979), the Hellas Planitia gravity low 
(Sjogren and Wimberly, 1981), and the Tharsis region (Janie and Erkul, 1991). 

Mars will receive its closest scrutiny starting in October of 1993, when the Mars Observer 
spacecraft arrives. Mars Observer has both a high-resolution camera (2.8 km feature resolution) 
and altimeter (30-m resolution) for surface characterization, and radio science gravity campaigns 
will be performed throughout the approximately 2-year mapping cycle. It is expected that from a 
50th to perhaps a 70th degree and order global field will be produced. Because of the nearly 
circular, low altitude (378 km) orbit and polar inclination of the spacecraft, the Mars Observer 
radio science team expects that the Mars gravity field will eventually be known more accurately 
than even that of the Earth. Discussions of gravity results from previous Mars missions as well as 
predicted Mars Observer contributions to Mars gravity field knowledge appear in Rosborough 
and Lemoine, 1991, and Tyler et al., 1992. 

Until the arrival of the Galileo spacecraft in 1995, Jupiter will have had only spacecraft flybys for 
gravity field investigations. Pioneers 10 and 1 1 both passed by the planet, closer in fact than the 
later Voyager 1 and 2 encounters. Data from all four flybys were combined in one study 
(Campbell and Synnott, 1985), and allowed mass calculations for the four Galilean satellites that 
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were a factor of 12 better than the analysis of Pioneer data alone. Even the combination of all the 
data was not sufficient to produce any harmonic terms for the satellite fields. Jupiter mass was 
determined to 0.00008%, and estimates of several of the low even harmonic coefficients were 
obtained (see Anderson et al., 1992c, for review). No features of the inner core of Jupiter or any 
other gas giant planet have ever been observed. 

The Galileo Orbiter is not expected to significantly increase understanding of the Jupiter gravity 
field, as its orbit will remain at a much higher altitude than that of any of the flyby missions. 
However, a great deal of new knowledge will be obtained regarding the larger satellites and the 
rotational and tidal characteristics of their interaction. For example, encounters with Io and 
Ganymede are expected to provide enough information to make decisions possible among 
different chemical and physical composition models (Anderson et al., 1992c). 

Both Voyagers (1980 and 1981), along with Pioneer 11, also made flybys of Saturn. Though 
Pioneer 1 1 flew closer to the planet, proximity of the ray path to the Sun introduced a large 
amount of plasma noise in the signal, and Voyager data were therefore an improvement (Null et 
al., 1981 for Pioneer results). The Voyagers also obtained closer looks at some of the larger 
Saturnian moons. A study combining the Voyager and Pioneer data was produced by Campbell 
and Anderson (1989), which also contains descriptions of previous results for Saturn gravitational 
characteristics. In it are improved masses for the moons Rhea, Titan, Tethys, and Iapetus, as well 
as several harmonic coefficients of Saturn's gravity field. In addition, upper limits on second-order 
coefficients of Titan's gravity field were also calculated. 

One small moon, 1981 S13, was first detected with the Voyager 2 Network-generated Doppler 
data. After gravitational perturbations were noticed in the spacecraft's motion, optical images of 
candidate body locations revealed the new satellite (Showalter, 1991). 

Saturn will be studied by the Cassini orbiter and the Huygens Titan Probe. It is hoped that 
Saturn's gravity field will be determined to a 2nd degree and order with possible higher zonal 
harmonics. As Cassini loops through the Saturnian system, encounters with several of the other 
large moons will fill in information on their masses and mean densities and better characterize the 
interdynamics of the planet and its satellites (Kliore et al., 1992). 

The best information on Uranus comes from the Voyager 2 flyby of the planet between late 1985 
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and early 1986. The combination of Doppler information and visual identification of mutual 
perturbations between the satellites allowed greatly improved mass determinations for the five 
major satellites, and therefore produced mean densities and indications of the moons' composition. 
The compositional information, in fact, led to the conclusion that the satellites could not be 
captured comets, as had been suggested (see Anderson et al., 1987a, for discussion). Uranus itself 
gained a more accurate value for both its own mass and the solar-to-Uranian-system-mass ratio, 
though significant information for other gravitational field parameters was not obtained. A model 
of Uranus' gravity field would be particularly interesting as Uranus is the only planet whose axis 
of rotation appears nearly parallel to the ecliptic plane. 

Neptune was the last encounter for Voyager 2, in mid 1989. Because of the close approach and 
the past difficulties in obtaining information about such a distant planet using only earthbound 
techniques, the Voyager flyby was able to increase the accuracy of the solar-to-Neptunian-system- 
mass ratio by a factor of 1000, using both navigation and Doppler data (Tyler et al., 1989). From 
the data, the masses of Neptune and Triton were estimated to 0.0003% and 0.3%, respectively. 
Even the small satellites Nereid and 1989N1 were observable in the spacecraft Doppler signature 
(Tyler et al., 1989). Estimates of a few, low order coefficients of the Neptunian gravity field were 
also obtained. A discussion of the interior structure ofNeptune appears in Hubbard et al., 1991. 

Pluto has yet to be approached by spacecraft, although there are studies for a possible future 
mission. As with Neptune, the increase in volume and accuracy of information from a spacecraft 
encounter over current Earth-based observations would be enormous. 

Most asteroids are too small to produce noticeable effects on spacecraft. Galileo's 1991 encounter 
with the asteroid Gaspra produced only images. In 1993, however, Galileo is scheduled to fly by 
the asteroid Ida, and it is expected that the relatively close approach to Ida will yield the first 
experimentally measured estimate of an asteroid's mass and therefore density. The corresponding 
constraints on its composition will allow new insight into theories of the evolution of the solar 
system and asteroid belt. 

It is interesting to note that over an 18-year period (1962-1980), the basic ground system 
Network-generated Doppler accuracy has improved 50 fold, from 5 mm/sec in 1962 to 0.1 
mm/sec in 1980. In 1992, Mars Observer is predicting a 0.01 mm/sec accuracy in overall ground 
system performance. 
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Figure 10. Doppler residuals and line-of-sight acceleration: (a) and (b) over Isidis Planitia, (c) and (d) over Olympus 
Mons; 1 Hz = 65 mm/second (from Sjogren 1979) 
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VI. Search for Low-Frequency Gravitational Radiation 


Gravity is described theoretically as a curvature of space-time. When a massive physical system 
undergoes violent dynamics, the curvature varies with time, radiating waves of space-time 
curvature (gravitational waves). Detection and subsequent detailed study of gravitational waves 
have the potential to revolutionize our view of the Universe. As emphasized by Thome ( 1987 , 
1 99 1 ), gravitational waves offer a "new window" for observational astronomy, giving information 
fundamentally different from that available in the electromagnetic window. The experimental 
difficulty with direct detection is that gravitational waves interact so weakly with matter. This 
extremely weak interaction will be an advantage in ultimate detections since the waves received 
will have propagated essentially unchanged from their sources; information about the time 
evolution of sources during violent events will have been preserved, unconfused by subsequent 
absorption or scattering. Relativistic motion of bulk matter and strong gravitational fields are 
central to theoretical views about violent activity in, for example, galactic nuclei and quasars. 
Gravitational waves from these objects will give us the first observations of the interiors of these 
strong-gravity, high-velocity regions. 

A gravitational wave is a propagating, polarized gravitational field-a "ripple" in the curvature of 
space-time. These waves are predicted by all relativistic theories of gravity. However, different 
theories of gravity predict different properties of the radiation (e g., propagation speed, 
generation mechanism, and polarization properties). The waves are characterized by their effects 
on separated test masses and separated clocks. The waves are propagating strain of spacetime, 
hence their amplitudes are dimensionless strain amplitudes (gravitational potential/c 2 ) that change 
the fractional difference in distance between test masses and the rates at which separated clocks 
keep time. 

In Einstein's theory of general relativity, gravitational waves are propagating solutions of the 
Einstein field equations. These waves carry energy and momentum. As with electromagnetic 
waves, gravitational waves are transverse to the propagation direction, have two independent 
polarizations, and propagate at the speed of light. However, unlike electromagnetic waves, 
gravitational waves couple to all matter (not just charged matter) and are extremely weak: the 
ratio of gravitational forces to electrical forces is about IQ- 40 . This extreme weakness has two 
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consequences: (1) gravitational waves are generated at detectable levels only by very massive 
sources under violent dynamical conditions (that is, astrophysical sources), and (2) the waves thus 
produced propagate unchanged from their sources, unaffected by scattering or absorption due to 
intervening matter. As mentioned above, this means that detections give information on the 
interiors of high gravity, high velocity regions not obtainable by other means. 


Alternate relativistic theories of gravity also agree on the existence of gravitational waves. Those 
theories differ, however, in essential features of the radiation: the number of polarization states, 
propagation speed, and efficiency of wave generation. 

Observationally, the strength of a gravitational wave is characterized by the fractional distance 
change between two test masses that is induced by the wave: the dimensionless strain 
amplitude," h=A^/f 0 , where i 0 is the unperturbed distance between the masses. 

Gravitational waves are conventionally classified on the basis of their temporal behavior, bursts 
are "on" for a well-defined time that is short compared to the observing interval, periodic waves 
have power concentrated at a few Fourier harmonics (having amplitude and phase that are 
sensibly constant over the observing interval), and stochastic waves have random variation over a 
time that is long compared with the observing interval. Hybrid waveforms can also exist, e g., 
quasi-periodic chirp waves. 

Bursts (which last for at most a few cycles) could be produced on a variety of time scales by 
collisions of stars or black holes, supernovae collapse to a neutron star, collapse of a star or star 
cluster to form a black hole, final coalescence of compact binaries (neutron stars and black holes), 
and the fall of stars or small black holes into supermassive black holes. 

Periodic waves (superpositions of one or more sinusoids that are approximately constant in 
amplitude and frequency over a typical observing time) could be produced by nonaxisymmetric 
rotating neutron stars, binary stars, and binary black holes. 

Stochastic waves (random fluctuations that persist for times long with respect to the observing 
interval) will be produced by an ensemble of incoherently radiating binary stars, deaths of 
pre-galactic massive stars, vibrations of cosmic strings, or, finally, as a relic of mhomogeneity in 
the structure of the big bang itself. 
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These waves are produced on a variety of time scales, depending on the masses of the objects 
involved. In general, stronger gravitational waves are expected to be produced on longer time 
scales. The time scales and strengths of the gravitational radiation from these sources can be 
estimated, based on theory and what is known about the universe from electromagnetic 
observations. Figure 12, adapted in highly simplified form from Thome (1987), shows the 
estimated strength and time scale of the radiation for some of these sources. It should be 
understood that these estimates are highly uncertain; the strength of the source's waves for a given 
distance from the Earth and the rate of occurrence of that type of source (thus the distance to the 
nearest one) are often extremely uncertain. This diagram should thus be viewed as current 
theoretical thinking, with the understanding that actual gravitational wave detections may provide 
surprises. See Thorne (1987) for a thorough discussion. 

Figure 12 shows the upper limit to the anticipated strength of burst sources based on current 
theoretical understanding. This limit is based on general physical considerations, including the 
assumptions that we do not live in a special time or place in the universe, that there are no 
enormous primordial bursts, that no single radiating object in the galaxy has a mass greater than 
100 million solar masses, and that there is no significant beaming of the radiation from the 
strongest sources. 

The "closure density" shown in figure 12 is the envelope of the family of curves having a 
bandwidth equal to a center frequency that represents the level of a stochastic gravitational wave 
background having energy-density sufficient to gravitationally close the universe. "Collapse to 
black hole" is the anticipated strength for burst radiation from black hole formation, assuming the 
mass of the hole (indicated) and assuming the source is at the Hubble distance (3 Gpc) or the 
Virgo cluster (where the amplitudes would be about 200 times larger). The shaded area labeled 
black hole binaries (Hubble distance)" indicates the amplitudes of periodic radiation from black 
hole binaries of indicated mass at the Hubble distance. The nearly horizontal arrows indicate 
evolutionary tracks of these coalescences at the Hubble distance. These would produce periodic 
and chirp radiation culminating in bursts at the time of coalescence. Amplitudes again scale 
upwards by about 200 for sources in the Virgo cluster. "Known galactic binary stars" shows the 
calculated strength of periodic waves from known binary stars in our galaxy. "Compact binaries" 
shows the calculated amplitudes for compact binaries, e g., neutron star binaries, in our galaxy. 
Loosely, this diagram resolves into three spectral bands: high frequencies (>10 Hz), low 
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frequencies (~10 -4 to 10 -1 Hz), and ultra-low frequencies (<10 7 Hz). Detection approaches are 
different in each band. 

Noise (notably seismic noise) causes ground-based detectors to be most sensitive to 
short-duration waves. There are two ground-based methods: resonant bars and laser 
interferometers (Michelson et al„ 1987; Abramovici et al„ 1992). Bar technology was pioneered 
in the 1960s by J. Weber of the University of Maryland. In a resonant bar antenna, gravitational 
waves excite mechanical oscillations in the bar; these oscillations are read out by a 
transducer/amplifier system. Current generation cryogenically cooled massive bars now achieve 
sensitivities of ~10 -18 . In laser interferometer antennas, the detector typically consists of three 
suspended masses arranged at the corners of a right-angle "L." Mirrors are attached to the masses 
on the ends of the arms and to a beam splitter on the corner mass. A gravitational wave incident 
on the system pushes the masses in one arm together while pushing the masses on the other arm 
apart. This is reversed on the next half cycle of the wave. The relative motions caused by the 
wave are detected as a fringe shift in a laser interferometer system monitoring the masses. The 
system is inherently broadband; many Fourier components are measured and thus the gravitational 
waveform can be constructed. The current generation of ground-based beam detectors has 
sensitivity comparable to that of resonant bars for -kHz frequency waves. The Caltech/MIT 
LIGO (Laser Interferometer Gravitational-Wave Observatory) collaboration is moving forward 
with construction of a pair of very large beam detectors (4 km arms) with more than three orders 
of magnitude better sensitivity expected. 

Isolation of the test masses from seismic noise, other acoustic noise from the environment, and 
fluctuating gravity gradients is prohibitively difficult. Therefore, for observations in the 
"low-frequency" (less than about 10 Hz), spectral band, the detectors must be space based. 

Space-based detectors are sensitive at low frequencies because of the large characteristic scale of 
the experiments (mass separations are approximately in astronomical units) and because of the 
low "environmental noise" in space. The current generation of space-based detectors uses the 
Earth and a distant spacecraft as free test masses. The Doppler tracking system of the Deep Space 
Network, driven by an ultra-high-quality frequency standard on the ground, monitors a 
transponded microwave link with a distant spacecraft. In doing so, it continuously measures the 
relative dimensionless velocity (Av/c) between the Earth and spacecraft. A gravitational wave 
passing through the solar system produces a strain of Af/^ 0 on the Earth-spacecraft system and 
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proportionately shifts the frequency of the clock driving the Doppler system. The result (figure 
13) is that the gravitational waveform is replicated three times in the Doppler tracking time series! 
once when the wave "buffets" the Earth, causing a small change in the difference between the 
transmitted and received Doppler frequencies; once when the spacecraft is buffeted by the wave, 
causing the signal transponded at an intermediate time to be different than the transmitted signal; 
and once when the initial Earth perturbation is transponded back to the Earth at a "two-way light- 
time" after the initial pulse. This three-pulse response is an important signature of a gravitational 
wave and allows discrimination against noise sources in the Doppler measurement system that 
have different response functions. 

The amplitudes of the three pulses are such that the areas of the three pulses add to zero. This 
defines the low frequency band edge in spacecraft observations: as the pulse duration becomes 
comparable to the two-way light-time between the Earth and the spacecraft, the pulses begin to 
overlap and cancel to first order. The high-frequency band edge is ~0. 1 Hz, set by a combination 
of thermal noise in the receiver and frequency standard and (for some spacecraft) high-frequency 
unmodeled motion of the spacecraft. 

At ultralow frequencies, pulsar timing can be used to search for stochastic waves from the early 
universe. The idea is that the rotation of the pulsar is a highly accurate "clock." If a gravitational 
wave is incident on the pulsar (or Earth), the rate at which the pulsar "clock" (or Earth clocks) 
keeps time is changed. Thus pulsar timing residuals can be used to constrain the spectral density 
of random gravitational wave perturbations at ultralow frequencies. For a review, see Section 5 
of Thome (1987), and references therein. 


The technique of relevance to the Deep Space Network is the Doppler technique. In this section, 
some of the details of this technique are given, and a discussion of the ultimate sensitivity of the 
technique is given. As outlined previously, the Doppler tracking measures the relative 
dimensionless velocity Av/c = Af/f 0 of Earth and spacecraft as a function of time, where Af is the 
perturbation in the Doppler frequency and f 0 is the nominal radio frequency of the link. A 
gravitational wave of strain amplitude h incident on the system causes small perturbations in the 
tracking record. These perturbations are of order h in Af/f 0 and are replicated three times in the 
Doppler data (Estabrook and Wahlquist 1975). That is, the "signal part" of the observed Doppler 
time series is the convolution of the waveform 
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g(t) = (1 - H 2 ) -1 {n • [h+(t) e + + h x (t) e x ] ■ n} 

with the function 

[(H - l)/2] 5 (t) - p 5[t - (1 + p)L/c] + [(1 + p)/2] 5 (t - 2L/c) 

where p is the cosine of the angle between the Earth-spacecraft vector and the gravity wave 
vector and L is the Earth-spacecraft distance, n is a unit vector from the Earth to the spacecraft, 
and e+, e x are transverse, traceless polarization tensors (Estabrook and Wahlquist 1975, 
Wahlquist 1987). The sum of the Doppler perturbations of the three pulses is zero; pulses having 
duration longer than about the one-way light-time produce overlapping responses in the tracking 
record, and the observed response then cancels to first order. The tracking system thus has a 
natural passband to gravitational excitation: the low-frequency band edge is set by pulse 
cancellation to be approximately 1/one-way light-time, while thermal noise in the radio system and 
short term stability of the frequency standard limits the high-frequency response to an order of 0.1 
Hz. 

Gravitational waves must be detected in the presence of noise. Experience to date indicates that 
the Doppler time series is adequately modeled as 

y(t) = Af/f 0 = gravity waves + propagation noise + clock noise + thermal noise + systematic errors 
= {[(p-l)/2]5 (t)-p5[t-( l+p)L/c]+[( 1 +p)/2)8 (t-2L/c)} * g(t) 

+ {5(t)+8[t-2(L)/c]} * troposphere(t) 

+ {8(t)+5[t-2(L-x)/c] } * plasma(t) 

+ {5(t)-5(t-2L/c)> * clock(t) 

+ thermal(t) 

+ systematic errors 

where * indicates convolution, g(t) is the scalar gravitational waveform (produced by the tensor 
wave's interaction with the Earth-spacecraft system), and "troposphere," "plasma," "clock" and 
"thermal" are noise processes, L is the distance to the spacecraft, and x is the distance to the 
equivalent solar wind "phase screen" (x * 0.25 AU for experiments done near solar opposition). 
Slowly varying systematic errors are, at current and near-future sensitivity levels, significant 
mainly for burst wave searches when the radio tracking beam is at low elevation angles. 

Note that the plasma and clock noises enter in essentially orthogonal ways and are different from 
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the three-pulse response. This was noted by Estabrook and Wahlquist (1975) and can be exploited 
to improve signal-to-noise ratio (sometimes dramatically) over what one might expect. (For a 
discussion of noise sources and their transfer functions to the Doppler observable, see Wahlquist 
et al., 1977, or Armstrong, 1989). 

Searches for gravitational waves in the low-frequency band have used the Viking, Voyager, and 
Pioneer spacecraft to determine upper limits for gravitational wave strengths (Anderson, 1977; 
Armstrong et al., 1979; Hellings et al., 1981; Anderson et al., 1984; Anderson and Mashhoon, 
1985; Armstrong et al., 1987; Estabrook, 1988; Anderson et al., 1990). All such searches have 
been limited at low frequencies by plasma propagation noise. The best sensitivities to date are for 
periodic waves: 7xl0~ 15 (one sigma) (Anderson et al., 1992b). 

Near-future observations will use the Galileo, Ulysses, and Mars Observer spacecraft (Estabrook, 
1988; Bertotti et al., 1992). Especially noteworthy is a coincidence experiment involving Mars 
Observer, Galileo, and Ulysses in the Spring of 1993. Expected sensitivity for periodic waves in 
that experiment is ~10 ^ (set by the X-band plasma scintillation noise on Mars Observer). 

In the late 1990s the Cassini spacecraft, with a Ka-band radio link and ultrastable ground and 
spacecraft electronics, will be used for much higher sensitivity observations. An engineering study 
(Riley et al., 1990) indicates that sensitivity to periodic waves could be ~10 -16 for one month of 
tracking near solar opposition, set by the ability to monitor and remove phase scintillations due to 
the wet troposphere. 

Sensitivity improvements in the low-frequency band significantly better than those possible with 
Ka-band and precision tropospheric monitoring will probably require moving both the test masses 
into space. This would remove the problem of tropospheric correction and the host of unmodeled 
geophysical effects (e.g., microseismic activity and unmodeled Earth tide effects) that become 
important at these very sensitive levels. One possibility, promising but not yet studied in detail, is 
to use a distant spacecraft coupled with a deep space station in a suitable Earth orbit as the two 
test masses. 
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Figure 12. Estimated strengths and gravitational wave radiation sources 



Gravitational 



Figure 13. Gravitational wave form three-pulse response in Doppler 








VII. Measurements of Relativistic Time-Delay 


One of the most exciting aspects of the capability to obtain range measurements from spacecraft 
in solar conjunction phases is a test of a portion of Einstein's theory of general relativity that 
predicts that electromagnetic waves passing close to massive bodies (e.g., the Sun) are delayed by 
the gravitational field. The signal time delay is dependent upon the observer-Sun-spacecraft 
geometry as follows: 


At = 2GM 0 c 3 (l+7) ln((r e + r m + R)(r e + r m - R) -1 ) 

where: 


At = signal delay, sec 

G = gravitational constant, km 3 /kg/sec 2 

M 0 = solar mass, kg 

c = speed of light, km/sec 

r e = Earth-Sun distance, km 

r m = Sun-spacecraft distance, km 

R = Earth-spacecraft distance, km 

In Einstein's theory of general relativity, the (parameterized post-Newtonian) parameter y is 
replaced by 1 Hence the nature of the experiment is to verify that y indeed equals 1 . 

To measure signal time delay, the Deep Space Network uses a binary coded sequential-acquisition 
ranging system. The hardware is called the Sequential Ranging Assembly (the predecessor to that 
was the Planetary Ranging Assembly). Ranging data are described in Section II. 

The first major attempt to measure the relativistic time-delay of a spacecraft signal occurred 
during the solar conjunction of Mariners 6 and 7 during mid-1970 (Anderson et al., 1975). This 
experiment was able to verify the value ofy to about 6%, i.e., y — 1.00 ± 0.06. The major error 
source that prevented that experiment from achieving greater accuracy was the non-gravitational 
forces on the spacecraft, which limited the accuracy of the orbit determination. Additionally, 
uncertainties in the signal delay caused by the scattering of free electrons in the solar corona 
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contributed an error of about 2% to the determination of y from Mariners 6 and 7. Later, non- 
gravitational acceleration noise was reduced by the Mariner 9 Mars Orbiter, and a measurement 
ofy to the 2% level was achieved (Anderson et al., 1978; Reasenberg and Shapiro, 1977). 

During 1976-1977, the relativistic experiment was performed using both the Viking Mars Orbiter 
and Lander spacecraft (Shapiro et al., 1977; Cain et al., 1978; Reasenberg et al., 1979). Three 
factors led to a very significant improvement in accuracy over the earlier Mariner experiments: 

(1) Significant improvement was made in the ground equipment; in particular the ability to 
calibrate the ground system ranging delay. 

(2) The availability of a landed spacecraft. The highly precise knowledge of the landed spacecraft's 
location eliminated the much larger positional errors of heliocentric spacecraft, such as those due 
to orbital errors and unmodeled forces. 

(3) Dual-frequency downlink (S- and X-band) capability. The dual-frequency downlink of the 
orbital spacecraft permitted accurate measurements of the signal delay due to electron density, 
thus significantly reducing the plasma uncertainty, which had previously been the largest error 
source. 

Independent determination of the parameter y at the Jet Propulsion Laboratory and Massachusetts 
Institute of Technology led to a value of 1 ± 0.002, the most accurate test of the radar time-delay 
to date (see Anderson et al., 1986, for a more detailed review). 

Since the Viking experiment, a test of the solar relativistic time-delay was performed with the 
Voyager 2 spacecraft during solar conjunction in December of 1985, shortly before the Uranus 
encounter in January of 1986 (Krisher et al., 1991b). The ranging system was similar to that used 
for Viking, with an S-band uplink and dual S- and X-band downlinks. For Viking, dual-band 
ranging was possible only with the orbiters, not the landers. Therefore, calibration of plasma delay 
to the landers could be done only indirectly by using the orbiter ranging data. With Voyager, it 
was possible for the first time to calibrate directly the plasma delay on the downlink. In addition, 
the conjunction was close to the solar equator, which minimized latitudinal variations in the solar 
plasma density (Anderson et al., 1987b). A remaining limitation in the Voyager experiment was 
non-gravitational accelerations of the spacecraft due to attitude control activity. This was 
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modeled by using telemetry records of thruster pulse firing. The final result was a determination of 
y = 1 ±0.03. 

Another relativistic time-delay test was attempted during the solar conjunction of Voyager 2 in 
1988. Because of the increased range to the spacecraft, three-way ranging was necessary for 
many passes. This provided an opportunity to calibrate accurately the plasma delay on the single 
frequency uplink by using the preceding, overlapping dual-band downlink measurements. 
Unfortunately, three-way ranging was still under development in the Network and much of the 
ranging data were corrupted by code-synchronization errors. These errors could be canceled by 
differencing the dual-band range measurements, thereby providing useful measurements of 
coronal electron density during solar maximum activity in the 11-year cycle (Krisher et al., 
1991a). It was not possible, however, to obtain an interesting test of the relativistic time-delay. 
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VIII. Spacecraft Deceleration in Cometary Comae 


The Giotto spacecraft encountered comet P/Halley on 13 March 1986 (Edenhofer et al., 1987) 
and comet P/Grigg-Skjellerup on 10 July 1992 (Patzold et al., 1993). The scientific objective of 
the Giotto radio science experiment was the determination of the cometary mass flux (gas and 
dust) on the spacecraft by observing its change in velocity during the flight through the cometary 
atmosphere (coma). In principle, the change in velocity due to cometary dust drag can be 
determined using three different radio navigational techniques. Two of these, based on two-way 
ranging and Doppler observations (S-band uplink, S- or X-band downlink), were conducted after 
the last pre-encounter spacecraft maneuver (two days before closest approach to the nucleus) and 
before the first post-encounter maneuver (one day after closest approach). During the actual 

encounter period, the X-band one-way downlink frequency (third technique) was recorded in real 
time. 


Comparison of the pre- and post-encounter two-way coherent Doppler and ranging residuals 
yields the total change in the spacecraft's radial velocity at encounter. These data were recorded at 
rates of one Doppler sample per second and one ranging point per ten minutes. Knowing the flyby 
geometry, one can compute the total change in spacecraft velocity along the relative flyby velocity 
vector and, from this result, the total impacted effective mass. Table 2 summarizes the results 
from the Halley and Grigg-Skjellerup flybys. 

The change in velocity as a function of time can also be extracted from the one-way noncoherent 
downlink frequency recorded during the actual encounter (closed-loop: 10 samples per second; 
open-loop: 50,000 samples per second). During the Halley flyby, it was obvious that the 
deceleration was not continuous, occurring rather in distinct steps interpreted as separate 
traversals of successive dust jets emanating from the nucleus surface. Unfortunately, a similar 
observation was not possible during the Grigg-Skjellerup flyby because the deceleration was 
below the noise level. Electrical and/or mechanical disruptions of the radio system were incurred. 
These were presumably caused by dust particle impacts on the spacecraft that generated 
frequency disturbances superimposed on the true Doppler frequency shift. Much larger Doppler 
shifts were observed in the one-way data than in the two-way data (16.9 Hz instead of 4.6 Hz at 
Halley; 124 Hz instead of 1 1 mHz at Grigg-Skjellerup). 
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Other changes in Giotto's dynamical state (nutation and spin rate) could be detected in addition to 
the deceleration (see Table 2). These were caused by impacts of very large dust particles (larger 
than 10 mg). The mass of the largest particle at Halley (about 400 mg) was estimated from a 
statistical analysis. At Grigg-Skjellerup, the nutation and spin rate changes could be most 
plausibly explained by one single dust particle of at least 30 mg out of the total impacted mass of 

39 mg. 
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Table 2. Comparison between Giotto flyby parameters at 
Comets Halley and Grigg-Skjellerup (GS) 


Parameter 

Hallev 

GS 

Relative flyby velocity, km/sec 

68.37 

13.99 

Angle between relative velocity 

and spacecraft spin axis, deg 

0 

68.8 

Angle between relative flyby velocity 

and direction to Earth, deg 

44.2 

67.3 

Distance to Earth, AU 

0.9 

1.4 

Distance of closest approach, km 

596 

<200 

Spacecraft mass, kg 

574 

529 

Change in radial velocity, mm/sec 

-167 

-0.4 

Total change in velocity, mm/sec 

-230.5 

-1.0 

Impacted effective mass, mg 

1932 

39 

Change in spin period, msec 

+12 

-1 

Nutation angle, deg 

1.04 

0.10 

Largest impacted particle, mg 

*400 

*30 
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